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Static and dynamical properties of elastic phase transitions under the influence of short-range 
defects, which locally increase the transition temperature, are investigated. Our approach is based 
on a Ginzburg-Landau theory for three-dimensional crystals with one-, two- or three-dimensional 
soft sectors, respectively. Systems with a finite concentration tid of quenched, randomly placed 
defects display a phase transition at a temperature T c (jid), which can be considerably above the 
transition temperature Tj? of the pure system. The phonon correlation function is calculated in 
single-site approximation. For T > T c (no) a dynamical central peak appears; upon approaching 
T c (no), its height diverges and its width vanishes. Using an appropriate self-consistent method, 
we calculate the spatially inhomogeneous order parameter, the free energy and the specific heat, 
as well as the dynamical correlation function in the ordered phase. The dynamical central peak 
disappears again as the temperatur is lowered below T c (nu)- The inhomogeneous order parameter 
causes a static central peak in the scattering cross section, with a finite k width depending on the 
orientation of the external wave vector k relative to the soft sector. The jump in the specific heat at 
the transition temperatur of the pure system is smeared out by the influence of the defects, leading 
to a distinct maximum instead. In addition, there emerges a tiny discontinuity of the specific heat 
at T c (no). We also discuss the range of validity of the mean-field approach, and provide a more 
realistic estimate for the transition temperature. 

PACS numbers: 64.60.-i, 61.72.-y, 63.20.Mt, 64.60.Ht 



I. INTRODUCTION 

The influence of defects on the statics and dynamics 
of structural phase transitions has been of considerable 
theoretical interest over the past two decades Jl]-^]. Es- 
pecially the appearence of a narrow central peak in the 
neutron scattering cross section, well above the transi- 
tion temperature, for both distortive Q and elastic struc- 
tural transitions ||l0[| , prompted various theoretical stud- 
ies dealing with local ordering phenomena around short- 
range static defects (for a review of the experimental 
facts, see Ref. and Ref. [O , and for a review of some 
theoretical results, see Ref. |13|). 

E.g., in Ref. 0] a one-dimensional model for continu- 
ous distortive structural transitions was studied, with the 
order parameter coupling to a single defect (Nd = !)• If 
the impurity locally increases the transition temperature 
T c ° of the pure system, this leads to a local condensa- 
tion of the order parameter in the defect vicinity. In 
higher dimensions, for this local order parameter conden- 
sation to occur, the defect potential strength must exceed 
a certain minimal threshold. Such locally ordered regions 
in the material emerging well above the pure transition 
temperature T® have played a prominent role in some 
of the theories attempting to explain the central peak 
phenomenon for distortive and elastic structural phase 
transitions H §. In this paper, we extend previous work 



on second-order ferroelastic phase transitions in d = 1 
H to higher space dimensions d, taking into account 
the crystalline anisotropy. To this end, we shall general- 
ize the methods developed for the distortive case ||] to 
(anisotropic) elastic systems, thus treating consistently a 
random impurity system with finite defect concentration 
jid (in the thermodynamic limit, both the number of lat- 
tice sites N — > oo and the number of defects Nd — > oo, 
but no = Nn/N = const.). 

In the framework of our mean-field approach, we shall 
find that defects which locally soften the crystal may in- 
duce a true phase transition at a temperature T c (nn) > 
T®. Below this defect-induced phase transition tempera- 
ture a spatially inhomogeneous order parameter emerges, 
whose average value remains very small in the vicinity of 
T c (nn) and only becomes noticeable near T®. Similarly, 
thermodynamic quantities (static susceptibility, specific 
heat, etc.) display prominent, broadened maxima near 
T® suggesting a "rounded" phase transition; however, the 
"true" singularities occur at T c (?id), but may not be seen 
in experiment at all, as their amplitude is only propor- 
tional to the defect concentration n^. The Bragg peaks 
of the low-temperature phase already appear in the scat- 
tering cross section for T < T c (nx>); as a consequence of 
the spatial inhomogeneity of the order parameter, they 
are accompanied by elastic Huang scattering peaks with 
finite q width. Furthermore, very close to T c (nu) an ad- 
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ditional dynamic central peak emerges, which may be in- 
terpreted as a dynamical precursor to the defect-induced 
phase transition. 

These mean-field results of course neglect order pa- 
rameter fluctuations, and exaggerate cooperative behav- 
ior. In reality, at T w T c (nr>) localized order parame- 
ter clusters appear, whose orientations however strongly 
fluctuate in space. Only at a lower temperature T or d < 
T c (ni)) (if at all) will they form a collective state with uni- 
form orientation, i.e.: the spatially inhomogeneous con- 
figuration predicted by mean-field theory. In order to 
provide a more realistic estimate of the proposed defect- 
induced transition temperature, we consider the cluster 
orientations as effectively Ising-like degrees of freedom, 
and then determine the cluster ordering temperature T or( j 
by calculating the free-energy difference of states with 
parallel and opposite orientation, respectively. Thus the 
onset of the order parameter, the Bragg peaks, and the 
Huang scattering will be shifted to somewhat lower tem- 
peratures, and the results of this work can essentially be 
used if T c (nu) is replaced by T or d. Provided that T or d 
is still considerably larger than Tj?, we thus expect the 
behavior of the thermodynamic quantities near T® to be 
very similar to the results presented here. 

This paper is organized as follows: In Sec. II we 
introduce the Ginzburg-Landau functional for a d— 
dimensional system, with one m-dimensional soft sec- 
tor, including randomly distributed point defects. The 
corresponding Langevin-type equation of motion is for- 
mulated. Furthermore we present an expression for the 
density-density correlation function, which serves as a 
starting point for subsequent considerations. In Sec. Ill 
the phonon response function is evaluated in the high- 
temperature phase, and the emergence of a dynamical 
central peak and a defect-induced phase transition well 
above T° is demonstrated. In Sec. IV we proceed to 
the ordered low-temperature phase, by using a suitable 
self-consistent approach. We determine the spatially in- 
homogeneous order parameter, the free energy and spe- 
cific heat, as well as the phonon correlation function, and 
discuss the singularities in these quantities. In addition, 
the scattering cross section S(k) (i.e., the density-density 
correlation function) is studied. In Sec. V we leave the 
realm of mean-field theory, and provide an estimate of 
the "true" defect-induced transition temperature (for the 
isotropic case), by identifying it with that temperature 
where already existing, but still fluctuating clusters con- 
dense to form a non-zero average order parameter. In 
Sec. VI, we briefly discuss the case of extended disorder 
(line or planar defects), and in Sec. VII we finally sum- 
marize and discuss our results. 



II. GENERAL EQUATIONS 



A. Model 

In order to describe elastic phase transitions of sec- 
ond order in d dimensions with an m-dimcnsional soft 
sector, we use an expansion of the elastic free energy of 
the unperturbed crystal with respect to phonon normal 
coordinates Q k , | |l4| , |l5[ . We disregard non-critical po- 
larizations; furthermore, aiming at the long-wavelength 
limit we keep only the lowest-order terms in the wavevec- 
tor expansion of the dispersion relation of the acoustic 
phonons. The wavevector k is then decomposed into its 
m-dimensional "soft" components p, and its (d — Tri- 
dimensional "stiff" part q, respectively: k = (p, q). Folk, 
Iro, and Schwabl have shown that terms of the form q 4 or 
q 2 p 2 are irrelevant (in the renormalization group sense) 
and do not affect the critical behavior of the system . 
In this spirit we use the following effective free energy 

F = J d d k J d d k' X - [(ap 2 + bq 2 + cp 4 ) S{k - k')Q k <2-k'] 

+0(Qi) . (2.1) 

The coefficient a is assumed to depend linearly on tem- 
perature, vanishing at T®: a = a'(T — T®); the very weak 
temperature dependence of the Ginzburg-Landau coeffi- 
cients b and c is neglected. 

In order to describe the influence of short-range de- 
fects, which locally increase the transition temperature, 
we assume that each defect creates a short-range poten- 
tial at its site, thus locally modifying the coefficients a 
and b of the Ginzburg-Landau functional ( |2.l| ) ; being in- 
terested in long-wavelength properties of the system, we 
can thus model the defect potential in the continuum by 
a S function. The coefficient a will be particularly sen- 
sitive to such a modification, as it becomes very small 
near the transition. For the coefficient c and the higher- 
order coefficients the defect influence is less important 
and will be neglected. We thus arrive at the following 
Ginzburg-Landau functional for the perturbed system, 

F = J d d k J d d k' X - [(ap 2 + bq 2 + cp 4 ) S(k - k')Qk<2-k' 

-0 k , k ,kk'Q k Q_ k ,] + 0(Qi) , (2.2) 

where </> k . k ' denotes the Fourier transform of the impurity 
potential (created by Nt> defects) 

JVd 

0(r)=C/ gj(r-r iD ) . (2.3) 

i D = l 

The defect strength U = a d X is taken to be positive, and 
therefore the transition temperature is locally increased 
at the impurities (here, V is the volume of the system, 
and Oq = V/N denotes the volume of the unit cell). 

The dynamics of the elastic crystal are governed by a 
Langevin-type equation of motion for the soft acoustic 
phonons Jig] , 
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Mco 2 Q k = - 



SF 



iMuj(Dp 2 + Dq 2 )Q k + r k + h k 



(2.4) 



The term on the left-hand side of Eq.(2.4) describes the 
acceleration, while the first term on the right-hand side 
provides the restoring force driving the system towards 
its equilibrium configuration. Note that we have intro- 
duced two different diffusive damping constants D and 
D for the soft and stiff sectors, respectively. r k denotes a 
stochastic force with vanishing average, (r k ) = 0; its sec- 
ond moment satisfies an Einstein relation, guaranteeing 
that exp(— F /k^T) is the equilibrium probability distri- 
bution. Finally, h is an external field which couples lin- 
early to the order parameter. Eq.(2.4) will be the basis 
for our discussion of the dynamical properties in the sub- 
sequent chapters. 



the following formula for the dynamical structure factor 
(for more details on the derivation, see Ref. ||) 



«S(k |W )» 



iv]>>, g + ]>>^<<^(k)}} 



a/3 



-2W 



2ttS(uj) + 



^fc Q /c /3 D Q ' 3 (k,t 

a/3 



-2W 



(2.8) 



The three different con tributions to the dynamical struc- 
ture factor in Eq.(2.S) are (i) the elastic Bragg peaks 
appearing at the reciprocal lattice vectors g of the actual 
crystal structure, given by the condition 



ig(a*+«^») 



1 



(2.9) 



B. Density— density correlation function 

In the following we shall primarily use a discrete lat- 
tice representation of the elastic system under consider- 
ation. The dynamic structure factor observed in scat- 
tering experiments is related to the Fourier-transformed 
density-density correlation function. Denoting the ther- 
modynamical average by (...), its definition is 



J \ 1 AS 



l<i,j<N 



(2.5) 



where denote the Bravais lattice sites (of the high- 
temperature phase), and the displacements from these 
equilibrium positions. 

In the discrete representation, with N lattice sites, we 
can write the Fourier-transformed defect potential as 



9kk< 



1 N ( Nn \ 

i,j=i \i D =i / 



-i(kxi-k'xj) _ (2.6) 



In a system with quenched, randomly distributed defects, 
all physical quantities have to be averaged over all possi- 
ble defect configurations [jl7). We denote this configura- 
tional average by ((■■■)); its formal definition reads 



JVd 

«...» =n 



i n 

-Y 



(2.7) 



In order to evaluate ((5(k, w))), we introduce a cumu- 
lant expansion for the combined thermal and configura- 
tional averages of e ik [ u ^°- ) ~ Ui ' t - | l and keep the terms up 
to second order. Next we decompose the deviations Uj(i) 
into a static contribution tpi and a fluctuating part Vi(t), 
and expand the exponential. Eventually one arrives at 



(ii) an additional static contribution to the structure fac- 
tor arising from elastic scattering from the random vari- 
ations of the local order parameter (Huang scattering) 

sf(k) = l£V*c»«-"i) - «V>?»«V#») ; 

(2.10) 

and (iii) the dynamical phonon-phonon correlation func- 
tion 



-ik{a i -a j )i,,at 1 \,3 



«?(*)<(0)> 



(2.11) 



which is connected with the dynamic phonon response 
function via the (classical) fluctuation-dissipation theo- 
rem 



D a/3 (k, 



2k B T 



Im G a/3 (k,w) 



(2.12) 



Finally, 

w = \ E kakP \(( w - (($*») (vf - «vf » 



a/3 

+ (((<*?) 



(2.13) 



is the Debye- Waller factor. Eq.(2.8) may be used for 
elastic as well as for distortive phase transitions. For an- 
tiferrodistortive transitions one has to sum over the dis- 
tinct sublattices in addition (see Ref. [pi) . The dynamic 
phonon-phonon correlation function Q2Tl| ) and the static 
Huang scattering contribution ( 2.1 0| ) will be discussed in 
more detail below. 
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III. HIGH-TEMPERATURE PHASE 

In order to calculate the phonon c orre lation function 
in the high temperature phase, Eq.( [2.2| ) is inserted in 
the equation of motion (2.4). Note that nonlinearities in 
the phonon normal coordinates are neglected, and thus 
fluctuations are only being accounted for in the Gaus- 
sian approximation. Upon differentiating the resulting 
expression with respect to hu 1 , and transcribing it to the 
corresponding discrete version, and finally using the fact 
that the average of the stochastic force vanishes, one 
arrives at the following mean-field recursion relation for 
the phonon response function 



Gkk' — GnirS) 



OkOkk' 



G O k^0kk"kk"G : 



k"k' 



(3.1) 



with the free phonon propagator 



G, 



Ok 



ap 



bq 2 



cp 



Id{m) = 



iMu(Dp 2 + Dq 2 ) . 

(3.2) 

(p 2 + q 2 ) 



Eq.( p.l[ ) may then be s yste matically iterated, and the 
configurational average ( [2.7| ) performed; e.g., applying a 
standard diagrammatic technique helps in collecting all 
the contributions to a certain order in the defect poten- 
tial strength A 0Jl^l ■ The configurational average yields 
a translationally invariant response function, 



((Gkk')) — Gk^kk' 



(3.3) 



Upon collecting all contributions which are proportional 
to the impurity concentration no = N^/N (single-site 
approximation; for more details, see Ref. ||), the result 
for the phonon susceptibility is 

G,; 1 = -Mlo 2 + ap 2 + bq 2 + cp 4 - iMuj(Dp 2 + Dq 2 ) 



Afc 2 n D 



1 - \(a /2Tr) d I d (m) ■ 
Here we have introduced the abbreviation 



o -Muj 2 + ap 2 + bq 2 + cp 4 - iMu(Dp 2 + Dq 2 ) 



d m p d d ~ 



Q 



(3.4) 



(3.5) 



m(d - m)T m Td- 



(p 2 +q 2 )f> 



-Mu 2 + ap 2 + bq 2 + cp 4 - iMuj(Dp 2 + Dq 2 ) 



dp dq 



and T n is the volume of the n-dimensional unit sphere 



T n/2 



r 3 



i 



(3.6) 



A denotes a natural short-wavelength cutoff (e.g., cor- 
responding to the Brillouin zone boundary), which also 
helps to ensure the convergence of the integral Id{m). 
Note that the dimensi on m of the soft sector in k space 
explicitly enters in Eq. (|3.5| ), and thus determines the im- 
portance of the fluctuation contributions. 

Eq.(3.4) implies the very remarkable result that due to 
the coupling to the softening defects, the entire system 
may become unstable towards a new ground state with 
finite average order parameter at a certain temperature 
T c (ni)), depending on the defect concentration jid; the 
criterion for this instability is 



lim [G k V = 0)/fc 2 







(3.7) 



As in the distortive case ||, we find that in general a 
certain minimal defect strength is required for this insta- 
bility to occur; yet, once this defect-induced phase tran- 
sition does exist, the associated transition temperature 
Tc(riY)) can be considerably higher than that of the pure 
system, Tj?. (In Sec. V, we shall comment on the validity 
of the mean-field approach, and estimate the transition 
temperature on a more realistic basis.) 

We have investigated three-dimensional systems with 
one one-, two- or three-dimensional soft sector, respec- 
tively. The qualitative features were found to be very 



similar in all these cases. The following figures refer to a 
three-dimensional system with a single one-dimensional 
soft sector. We have tried to use model parameters 
appropriate for NbjSn, which displays a second-order 
elastic phase transition near T — 45K [^(J; accord- 
ingly, we have used numerical values calculated from 
Refs. |0||§ (Table I). Thus, we have taken T c ° = 45K, 
no = 10 , • and adjusted the defect strength in order that 
T c (ni)) = 65K. However, a few remarks are in place here 
to explain some sources of inaccuracies. The assumption 
that the Ginzburg-Landau parameter a is merely linearly 
temperature-dependent is valid only near the phase tran- 
sition temperature of the pure system T®. Furthermore 
we approximated c and b as independent of temperature, 
and in addition assumed b to be independent of the direc- 
tion of the k vector in the stiff plane. This is not generally 
the case for NbsSn, but appears to describe the critical 
region well. The numerical values of the diffusion con- 
stant D — D and the nonlinearity d (see below) had to 
be estimated without reference to any experiment. 

Fi g.l d epicts the phonon correlation function D(k, ui) 
[Eq.( 2.12] )] for different temperatures T > T c (nD), evalu- 
ated for several angles 9 between the external wave vector 
and the soft sector. As becomes apparent in Fig.l, a dy- 
namical central peak in the phonon correlation function 
emerges in addition to the soft phonon peak (compare 
Ref. Q| for the one-dimensional case) . The height of the 
central peak grows, and its width decreases as T c (nn) is 
approached. The intensity of the central peak decreases 
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upon increasing the angle 9 between the wave vector k 
and the soft sector. This reflects the fact that wavevec- 
tors in the stiff sector do not probe the critical properties 
of the material. The dynamical central peak may thus be 
understood as a dynamic precursor to the defect-induced 
second-order phase transition at T c (riT)). 



IV. LOW-TEMPERATURE PHASE 

In this section, we use a self-consistent approach de- 
signed for the calculation of the order parameter v, the 
specific heat, the phonon correlation function and finally 
the dynamical structure factor in the ordered phase, i.e., 
for T < T c (nu)- 



A. Order parameter 

The starting point for the calculation of the order 
parameter is the full nonlinear Ginzburg-Landau func- 
tional, which in the discrete lattice representation reads 

N , JV JV D N 



N 

2=1 



i—1 *d = 1 



(4.1) 



Here denotes the value at lattice site i of that combi- 
nation of strain tensor components serving as the order 
parameter for the transition, hi is the corresponding ex- 
ternal stress acting on site i, and the static propagator 
Goij is defined by its Fourier transform 



G 0k 

n -i _ ap 2 



bq 2 + cp 4 



for d > m. 



(4.2) 



In the framework of the Ginzburg-Landau approxima- 
tion, i.e., neglecting order parameter fluctuations, the 
following stationarity condition can be derived (with 
hi = h = const.) 



i.e.: we assume that the order parameter at each lattice 
point i may be written as the sum of a homogeneous 
background A and an additional contribution B, if there 
is a defect at site i, thus enhancing the total value of the 
order parameter to A + B at the defect sites. Thus we ex- 
plicitly assume that at all defect sites the order parameter 
points in the same direction, and in addition neglect the 
spatial variation of the order parameter near the defects. 
However, as we shall see shortly, the second, seemingly 
very crude approximation already contains the possible 
relevant modifications caused by the impurities, namely 
(i) an enhancement of the spatially averaged order pa- 
rameter (corresponding to the parameter A), and (ii) the 
ensuing "screening" of the defect potential (described by 
the coefficient B). The more stringent approximation is 
the uniform orientation of the defect clusters, as implied 
by the mean-fiel d ap proach (see Sec. V). 



Inserting Eq. (4.4) into the stationarity equation (4.3) 
yields the recursion relation 



f k = /i<5 ko Go(k) + Go(k) y~] <?W^k 



k' 



(4.5) 



where we have introduced a renormalized propagator 

Go(k)- 1 =G (k)- 1 +dA 2 , (4.6) 

and a screened defe ct potential (f)kk' with weakened 
strength [see Eq.(2.6)] 



\ = \-d[{A + B) 2 -A 2 } 



(4.7) 



From Eq.(4.5) and the averaged stationarity equation 
we may derive two coupled nonlinear equations that 
uniquely determine the mean order parameter: (i) It- 
crating the recursion relation (4.5) in a similar way as 
for the dynamics in the previous paragraph, performing 
the configurational average, and summing the single-site 
contributions, one arrives at 



dA 2 - 



An D 



1 - \{a /2Tr) d J d 



(4.8) 



with the abbreviation 



SF 



= O G 0ij v i ~ A vA ^ 



dvf = h . 



(4.3) 

A general solution of Eq.( |4.3| ), with its combined non- 
linearity and randomness, poses a difficult problem. We 
thus use an additional approximation, namely the fol- 
lowing ansatz || for the thermodynamical average of the 
order parameter (denoted by P), 



A + Bj26 itin , 



(4.4) 



Jd 



d d k 



{a + dA 2 )p 2 + (b + dA 2 )q 2 + cp 4 



(4.9) 



(ii) On the other hand, immediate averaging of Eq. ( |f.5| ) 
yields 

(a + dA 2 )(A + n D B) - \n D (A + B) - h = . (4.10) 



Very assuringly, Eqs.(4.8) and ( |4.1C| ) yield non-zero solu- 
tions for ((£")) precisely below T c (ud) as determined from 
the high-temperature phase. Fig. 2 shows that the order 
parameter of the perturbed system as function of T looks 
similar to the corresponding curve for the pure system, 



5 



with the singularity at T® being smeared out by the disor- 
der. The order parameter sets in continuously at T c {nr>), 
with the usual mean field exponent (3 = 1/2, assumes 
small but finite values in the range T c (fiD) > T > T c °, 
and starts to grow to larger values only in the vicinity of 
T c °. Thus the transition temperature of the pure system 
remains an important parameter even in the perturbed 
system, while the true phase transition at T c (riD) may 
in fact be hardly noticeable in experiments. As before, 
the results for a three-dimensional system with one one- 
dimensional soft sector are depicted, but the qualitative 
features remain essentially the same in the cases of a two- 
or three-dimensional soft sector. 



In Fig. 5 the static phonon susceptibility G(k) = 
G k (a; = 0) (i.e., the inverse elastic constant, as modi- 
fied by the defects) is shown. The small but sharp peak 
at T c (n-£,) reflects the preordering of the defect regions, 
while the broad and much more prominent peak near 
yo corresponds to the ordering of the pure bulk crystal, 
though under the influence of the randomly spaced fields 
originating from the defect clusters; compare Figs. 2 and 
3. 



D. Dynamical structure factor 



B. Specific heat 

From the knowledge of the mean order parameter, we 
can rea dily calculate the (averaged) free energy from 
Eq.(l.l) in Landau approximation, and via 



G„ 



-T 



fd 2 F\ 
\dT 2 ) v 



(4.11) 



derive the specific heat C v , see Fig. 3. Obviously, the dis- 
continuity at T® has been smeared out, in place of which 
a tiny jump emerges at T c (n^). Although the phase tran- 
sition clearly occurs at T c (jin), the transition tempera- 
ture T® of the pure system remains of considerable im- 
portance; e.g., there is a distinct maximum of the specific 
heat near Xj?, while the extremely minute jump at T c (nu) 
might not be experimentally detectable at all. 



C. Phonon correlation function in the ordered phase 



In order to find the phonon correlation function in the 
temperature region with a finite order parameter, one 
again has to start from the full Ginzburg-Landau func- 
tional and use the ansatz for the order parameter (4.4). 



The crucial point is that one may then absorb the nonlin- 
ear term of the equation of motion in modified coefficients 
of the linear terms as follows 



a -> a + 3dA 2 



b^b + MA 2 , A -> A - UB{2A + B). 

(4.12) 



With these modifications one can use the same equations 
as in the high-temperature phase. 

The result is depicted in Fig. 4. The central peak in the 
correlation function disappears again when the temper- 
ature is lowered below T c (nT>). This dynamical central 
peak is thus confined to the region around T c (nn). As 
in the high-temperature phase, the intensity of the cen- 
tral peak decreases upon increasing the angle between 
the external k vector and the soft sector (with fixed tem- 
perature) . 



In order to describe scattering experiments, we have 
to calculate the density-density correlation function of 
Sec. 2. The first term in Eq.(2.8) yields the Bragg scatter- 



ing, and does not require any further comment; the third 
term is connected with the phonon correlation function, 
and has been discussed in the previous subsection. We 
therefore turn our attention to the second term. Tak- 
ing into account the soft acoustic phonon mode only, as 
above, we have t o cal culate the configurational average 
of fc 2 i/ , kV , -k [Eq.(2.10)]; using the same approximat ion s 
as in the beginning of Sec. 3, we may use equation ( |4.5| ) 
in the form 



fcV>k = frG (0)<5 k o + ^o(k) ^2 <^kk'^' 



(4.13) 



which yields 

k 2 M- k = /i 2 G (0) 2 5 ko + /iGo(O) 2 0ok'*>k'<5ko 



G (k)y"](/> kk /ivi/_ k 

k' 



(4.14) 



For this equation again a diagrammatic representation 
can be derived and in single-site approximation (i.e.: 
to order no) we find the following result (kx denotes the 
components of the wave vector which are parallel to the 
polarization of the soft mode) 



fc^(k) = n D \ 2 ((D}} 2 (k L ) 2 Go(k)G (-k) x 

-1-2 



(4.15) 



This expression can be further reduced using Eq. (^Ja). 
Finally, we arrive at 



fc 2 S c LL (k) = (a + dA 2 ) 



n D 



(k L ) 2 G (k) 2 . (4.16) 



Collecting all results, the final expression for the dy- 
namical structure factor reads 



G 



«S(k, W )» = 



(fc L ) 2 ((0) 2 /s 



A- 2 



"D 



Go(h) 2 (a + dA 2 Y 



-2W 



2ttS(uj) + (k L ) 2 D(k,oj)e 



-2W 



(4.17) 



Thus we have found three distinct effects, nameiy (i) 
new positions of the Bragg peaks as a result of the fi- 
nite order parameter [shifted and possibly new recipro- 
cal lattice vectors, see Eq.(2.9)] ; (ii) Huang scattering 



as a result of the spatially inhomogeneous order param- 
eter configuration, leading to a static central peak with 
finite width 7 = y/Ja~+ dA 2 )/c (in the soft sector) in 
Fourier space; and (iii) inelastic scattering, described by 
the phonon correlation function. Fig. 6 shows how the in- 
tensity of the Huang scattering varies with temperature 
for different wave vectors k = p in the soft sector. This 
additional elastic contribution sets in at T c (jid), and then 
grows to considerable values near Tj?. 



V. ESTIMATE OF THE CLUSTER ORDERING 
TEMPERATURE 

All our previous results for the statics were based en- 
tirely on the Ginzburg-Landau approximation, and dy- 
namic quantities were calculated in the Gaussian ensem- 
ble. This mean-field treatment of course neglects fluc- 
tuations, and apart from the fact that the critical ex- 
ponents will be changed near the transition, we have to 
consider the possibility that the above described defect- 
induced phase transition at T c (ni)) will disappear when 
fluctuations are properly taken into account. Namely, 
our mean-field approach basically implies that as soon 
as local condensates form near the defects, they immedi- 
ately lock into some cooperative state and form a non- 
vanishing average order parameter. In reality, proba- 
bly first these clusters may emerge at the defect posi- 
tions, however still quite independently fluctuating be- 
tween their different possible orientations. Only as the 
temperature is lowered even further, they will form a 
collective vibrational mode which finally condenses to a 
static order parameter at the "true" transition temper- 
ature T or d, with T c ° < T or d < T c (nn), the mean-field 
transition temperature. One would expect that such col- 
lective behavior of the distinct localized order parameter 
clusters arises when the correlation length of the pure sys- 
tem £, which determines the size of the defect-induced 
condensates, becomes of the order of the average defect 
separation rrj- A somewhat more favorable estimate re- 
sults from the argument that it should actually suffice 
when £ becomes large enough such that the distinct con- 
densates form a percolating cluster throughout the sam- 
ple; the condition for cooperative behavior then becomes 
£ « (n c ) 1 / d rD, where n c denotes the percolation thresh- 
old. 

In the following we give a more precise estimate T or( j , in 
order to see if it may still be considerably above the tran- 
sition temperature of the pure system T®. Our strategy 



is to calculate the free-energy difference AF between the 
following two configurations in a two-defect system be- 
low the temperature where localized clusters may form in 
two different orientations: (i) both order parameter con- 
densates oriented in the same direction, and (ii) oppo- 
site condensate orientations. The ensemble of localized 
clusters can then be effectively mapped onto an Ising 
system, with AF assuming the role of the exchange cou- 
pling. The critical temperature is now readily estimated 
as ksTord ~ AF/a,Q (gsq is the volume of the elementary 
cell) . We emphasize that we shall restrict ourselves to an 
isotropic system here, and consider the general case of 
an order parameter f described by the usual Ginzburg- 
Landau expansion of the free energy, which rather cor- 
responds to the case of distortive structural transitions, 
as studied in Ref. ||. However, the qualitative behavior 
is expected to be very similar for the anisotropic elastic 
phase transitions. 

Using the continuum representation, the free energy 
for a system in three dimensions with a single defect in 
the origin reads EJ 



F 



d 6 r [a-(b(r)]f(r) 2 +c[\7f(r)Y 



(5.1) 



where </>(r) is the positive 8 function defect potential with 
strength U — <ZqA. The stationarity equation then be- 
comes 



cVV(r) = [a - (/>(t)](p(t) + df(rf , 

which for a > may be approximately solved by 

,r/S p -r/£ 



(5.2) 

(5.3) 
(5.4) 



c/a is the correlation length of the pure 
(4ttR 3 X/3 - 10a) /d, and 
R 3 X = 1207rc 3 (C/- 1 - U' 1 ) 2 , with U m = 2irca denot- 
ing the minimum defect strength required for the local 
order parameter condensation to occur. 

Using these results, we can proceed towards the two- 
defect system with </>(r) = US(x)S(y)[5(z — Hd/2) + S(z + 
rrj/2)] by a simple linear superposition ansatz; i.e., we 
shall evaluate the free energy difference between the 
states 




system for T > T c °, f\ 



<P± = v{ x , V^ z - ? 'd/2) ± <p(x, y,z + r D /2) 



(5.5) 



By inserting into Eq. (|j) one readily finds the defect 
contribution 
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(5.6) 



as well as the linear overlap integral (conveniently eval- 
uated using elliptical coordinates) 



AFiin « Wnapl^e 



(5.7) 



the nonlinear overlap integral turns out to be of order 
e -2r D /£ anc j can thug b e neglected for £ < ro, when 
compared to the previous terms. Hence we find for the 
required free-energy difference 



AF « 167rc<^£e _rD/s 



I - ^ 

27rcrD / 



(5.8) 



using the above numerical values, we see that the defect 
contribution can in fact be neglected here. 

Hence we arrive at our final estimate for the cluster 
ordering temperature, which we identify with the "true" 
defect-induced transition temperature 



feTord ~ 167rav2 (£/a ) e 



(5.9) 



This expression may be cast into a somewhat more ex- 
plicit form by observing that the average defect separa- 

— 1/3 

tion can be written as td = aon B ; thus the required 
defect concentration for the transition to occur at a cer- 
tain value T = T or d becomes 



n D = — In 



167ra r 



knT, 



B^ord 



«0 



(5.10) 



from which T or d as function of tid may be inferred by 
inversion [£ c = \fcfa c 



a c = a'(T ord - T c )]. The result 
is depicted in Fig. 7. It can be seen that the calculated 
cluster ordering temperature may indeed be considerably 
above the phase transition temperature of the pure sys- 
tem T c ° = 45K, however, much larger impurity concen- 
trations riD are required than in the previous mean-field 
analysis. 

For the anisotropic elastic systems discussed in the 
bulk of this paper, fluctuations will be even less impor- 
tant. We conclude this section with the remark that the 
upper critical dimension as function of the dimension m 
of the soft sector was found to be @ 



deM = 2 + — ; 



(5.11) 



thus in three dimensions mean field theory yields exact 
results for a system with a one-dimensional soft sector, 
while for the case of a two-dimensional soft sector merely 
logarithmic corrections are to be expected. 



VI. EXTENDED DISORDER: LINE AND 
PLANAR DEFECTS 

We now return to the case of elastic phase transitions, 
and address the question of the influence of extended de- 
fects in contrast to the previously treated point disorder. 



Our system now contains randomly placed, but parallel 
linear or planar defects; the accor ding ly modified corre- 
lated defect potential (compare Eq. |2.3| ) , reads in the case 
of line disorder 



(r) = U ^ S(x - x lD )5(y - y io ) , 



(6.1) 



where x,y,z are the components of r, and z denotes the 
direction parallel to the lines; id labels the line de- 
fects. The defect potential for planar defects is defined 
analogously, namely for planes normal to the x direction 



A'd 



(r) = U ^ 8{x 



(6.2) 



With these definitions the same calculations as before 
may be performed, and it becomes obvious that the for- 
mer integrals reduce to integrals over the k vectors per- 
pendicular to the defects. One gets qualitatively the same 
results as in the case of point defects. 

In order to compare the effect of the different kinds 
of defects, we have calculated the order parameter in all 
three cases (points, lines, and planes) for the same de- 
fect strenght and the same defect concentration (i.e., the 
extended defects are viewed as correlated accumulations 
of point defects with the total number of - pointlike - 
defects held fixed). Therefore, the resulting differences 
solely originate in the different disorder dimensionality. 
The one-dimensional soft sector was taken to be perpen- 
dicular to the defects in order to provide a meaningful 
comparison. The result is depicted in Fig. ||; it can be 
seen that the effect of the defects is not a monotonous 
function of their dimensionality, but depends on the 
strength of two competing effects. On the one hand, 
when the temperature is lowered towards the phase tran- 
sition temperature and the correlation length £ grows 
accordingly, the order parameter cluster around a d'- 
dimensional defect grows proportional to . This 

effect renders low-dimensional defects more effective in 
influencing bulk properties. On the other hand, the sys- 
tem has a finite stiffness, characterized by the parame- 
ter c; and as a system with uncorrelated defects is more 
inhomogeneous than one with the identical amount of 
correlated disorder, this effect favors high-dimensional 
defects, because then the system stiffness can be more 
easily overcome by the joint action of neighboring de- 
fects. With the specific numerical values we have used, 
the line defects have only a tiny effect on the order pa- 
rameter curve in comparison with the point defects. The 
effect of the planar defects lies in between. We empha- 
size that this scenario could be different for other values 
of the stiffness parameter c. Finally, we remark that if 
one performs the above calculations for a system with 
extended defects, where the soft sector is not perpendic- 
ular to the defects, additional angle dependences ensue, 
and one has to add an additional term of the form cq 4 
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to the functional ( |2.2| ), in order to correctly account for 
the stiffness, which tends to prevent the building-up of 
order parameter clusters. 



VII. SUMMARY AND DISCUSSION 

In this paper we have studied the influence of point 
and extended defects on a d-dimensional elastic system 
with an m-dimensional soft sector undergoing an elastic 
phase transition of second order. We have calculated 
the phonon-phonon correlation function in the high- 
temperature phase. At a certain temperature T c (riD) 3> 
T®, an instability marking a defect-induced phase tran- 
sition may emerge, if the defect potentials are sufficiently 
strong. Above T c {n^) a dynamical central peak emerges 
in the phonon correlation function, whose intensity grows 
as the temperature is lowered towards T c (nrj). Contrary 
to the case of distortive phase transitions M, the max- 
imum of the central peak is exactly at u> = for all 
temperatures and not at small but finite frequencies . 
This, however, does not imply that the acoustic impu- 
rity modes are localized; at least for a simplifying one- 
dimensional single-defect model no additional localized 
impurity modes appear, but the defect rather causes 
a localized vibrational contribution to the propagating 
scattering states; for the long-wavelength phonons this 
quasi-resonant vibration then condenses at T c (jib) and 
forms the local order parameter clusters |7| . The dynam- 
ical central peak may be regarded as precursor of this 
phase transition; its height also depends on the angle 
between the external wave vector k and the soft sector. 
The smaller this angle, the more pronounced is the cen- 
tral peak. 

In the low-temperature phase T < T c (nr>), where a 
finite, and be it ever so small, order parameter exists, we 
have used a self-consistent mean-field calculation in or- 
der to calculate the average order parameter, the free en- 
ergy and specific heat, and the phonon correlation func- 
tion. The order parameter sets in continuously at T c (nn), 
remains very small in the temperature range between 
T c {nu) and T°, and reaches appreciable values only near 
T c °. In this way the order parameter curve resembles 
a somewhat rounded curve of the pure system. Analo- 
gously, the temperature dependence of the specific heat 
looks like the corresponding smeared-out curve for the 
pure system. The jump at T® is rounded, and a minute 
jump at T c (ni)) appears. Thus the phase transition of the 
perturbed system no longer occurs at Tj? but at T c (?7.d)- 
However, the phase transition temperature T® of the pure 
system remains important, as the remnants of the pure 
transitions induce marked, but rounded maxima in quan- 
tities like the order parameter susceptibility or the spe- 
cific heat near T c °. 

Having thus determined the order parameter, we were 
able to calculate the phonon correlation function in the 
ordered phase. The dynamical central peak disappears 



again as the temperature is lowered below T c (ni)). The 
dependence of the central peak on the angle between the 
momentum transfer vector k and the soft sector is very 
similar as above T c (71d)- The density-density correla- 
tion function determining the cross section for scatter- 
ing experiments consists of three terms: first, the term 
describing elastic Bragg scattering, second, a term cor- 
responding to Huang scattering caused by the spatially 
inhomogeneous order parameter configuration; this term 
yields a contribution to elastic scattering, leading to a 
static central peak with finite k width. The third term 
finally describes inelastic scattering and has been dis- 
cussed along with the phonon correlation function. 

We have also discussed the validity of our mean-field 
approach and provided an estimate (in the isotropic case) 
for the "true" precursor T c , defined as the cluster or- 
dering temperature T or d, which has to be distinguished 
from the temperature T c (tid) where localized, but fluctu- 
ating order parameter clusters appear. Only below T or( j, 
the previously independent ordered regions form a col- 
lective state leading to a non-zero average order param- 
eter. Generally, cooperative behavior of the defects is 
to be expected when the correlation length of the pure 
system £ becomes of the order of the mean defect separa- 
tion rt>. Although the ensuing cluster ordering temper- 
ature is considerably lower than T c (71d), the qualitative 
features of the present theory should remain largely un- 
affected, provided T c (tid) is replaced by T or d; i.e., the 
mean order parameter appears rounded near T c °, while 
static order parameter susceptibility and the specific heat 
display a strong but broadened maximum there. Fur- 
thermore, in the anisotropic systems under consideration 
here, fluctuations may actually be suppressed, rendering 
mean-field theory more reliable. Finally, the time scale 
of the order parameter condensate fluctuations will di- 
verge oc (T — Tord) -1 upon approaching T or d, which in 
experiment would eventually render them indistinguish- 
able from static inhomogeneities, leading to quasi-elastic 
Bragg and Huang scattering peaks. We have also inves- 
tigated the case of parallel extended defects (lines and 
planes), and found essentially the same features |2l| . 

These considerations led us to the qualitative phase 
diagram displayed in Fig. 7. For very tiny disorder con- 
centrations, the picture of isolated defects applies. Al- 
though preordered clusters may form considerably above 
T®, the ensuing condensates fluctuate independently and 
do not form a state with nonzero average order param- 
eter. The cluster reorientation rate will become very 
low as Tj? is approached, and coupling of these slow 
modes to the soft (acoustic) phonons will then lead to 
a dynamical central peak, see Ref. For higher, but 
still small defect concentrations, the clusters emerging 
at T c (n rrn D) will form a state with preordered defect 
regions and finite, but small average order parameter 
below T or d- This phase transition leads to discontinu- 
ities in thermodynamic quantities, like the specific heat 
or the static susceptibility, which are, however, probably 
unnoticeably small in experiment. On the other hand, 
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upon approaching T®, the phase transition temperature 
of the pure crystal, the bulk system orders; due to the 
influence of the preordered defect regions, the transition 
temperature will be slightly higher than Xj?, and the for- 
merly sharp singularities of the mean order parameter, 
specific heat, and static susceptibility appear character- 
istically rounded. The above described theory applies 
precisely to this concentration range (for very strong dis- 
order the single-site approximation breaks down). The 
central peak phenomenon is thus explained by a com- 
bination of clastic Bragg peaks of the low-temperature 
phase and static Huang scattering with finite width in q 
space. 

However, a different scenario is also conceivable, 
namely that as the cluster reorientation times become 
very long, the different, still independent condensates 
freeze in with spatially fluctuating orientations. The en- 
suing configuration would constitute a metastable state 
which is separated from the true ground state by high 
free energy barriers A. The typical flip rate would then 
be proportional to exp(— A/k^T); therefore, at low tem- 
peratures the true thermodynamic ground state may not 
be reached. In the spirit of the discussion in Sec.V one 
could possibly map this problem onto a random-field 
Ising model (for recent reviews, see, e.g. Ref. [0), the 
lower critical dimension of which is di = 2, and therefore 
long-range order is not destroyed in three dimensions. 

At last, we would like to contrast our picture with that 
of "glassy" systems. Although some of the features of the 
dynamical structure factor in glasses are at least quali- 
tatively similar, e.g., an elastic peak with finite q width 
appears, and the static susceptibility and specific heat 
may display characteristically rounded and broadened 
maxima, there are important differences. First, the or- 
der parameter of the pure crystal v would not constitute 
an appropriate order parameter for such a disordered, 
glassy system. Second, the character of the phase tran- 
sition should be entirely different, and in fact lead to ex- 
perimentally distinguishable behavior. The scenario de- 
scribed here is a genuine second-order phase transition, 
though induced by disorder (which locally softens the sys- 
tem); i.e.: critical phenomena are confined to regions very 
close to T c (tid) (or T or( j), and to wave vectors in the soft 
sector with p w 0. The freezing-in into a glassy state, 
on the other hand, would have to be described by an 
crgodicity-breaking (Edwards-Anderson) order parame- 
ter, and should actually be rather insensitive to k. Such a 
glass instability occurs, e.g., in orientational glasses |p3| , 
and possibly in relaxor ferroelectrics j24|]. We finally re- 
mark that for the case of first-order martensitic trans- 
formations, a model with disorder of the random T c type 
has been proposed, which can then be mapped onto a 
spin glass, and the ensuing glassy features were suggested 
to explain the prominent tweed microstructure found in 
these materials, as well as the central peak phenomenon 
there §|. 
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TABLE I. Ginzburg-Landau parameters, as used in the figures, if not specified otherwise [a = a'(T — T° )]. 

T° = 45 K 

a' = 1.0771 • 10- 12 erg K" 1 
d = 1.1363- 10 _8 erg 
M = 1.314 ■ 10~ 21 g 
b= 1.570- 10~ 10 erg 
c = 5 • 10~ 26 erg cm 2 
N/V = 6.751 ■ 10 21 cm- 3 
A = 3.0994 ■ 10" n erg 
D = D = 1 • lO^cmV 1 
k = C- V2-a* 

a* = 2-K/a = 1.189 • 10 s cm" 1 

T c (n D ) = 65K 

n D = 1 • 10~ 5 

A = 3.0994 ■ 10 _11 erg 
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FIG. 1. Phonon correlation function D(k,cj) = 2fcBTImG(k, w)/w [in 10 _27 cm 2 s] vs. phonon energy [in meV] for different 
temperatures and fixed angle 6 = (a), and for fixed temperature T = 65.01 K with different angles 9 (b); £ = kao/2 3 ^ 2 Ti = 0.02 
was used. 
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FIG. 2. Average order parameter {{v}} vs. reduced temperature t for iid = 10~ 5 , A = 3.0994 • 10 _11 erg (a), and «d = 10~ 3 , 
A = 3.0949 • 10~ n erg (b). The temperature range near t c {nu) = 0.444 is displayed in the insets. Note the different scales of 
the two insets; the disorder strength was adjusted to yield the same T c (nr>) in both cases (a) and (b). 
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FIG. 3. Specific heat C v vs. reduced temperature t for no = 1CP 5 (a) and no = 1CP 3 (b). The defect-induced temperature 
has been adjusted to t c (no) = 0.444 in both cases by changing the defect potential strength accordingly, see Fig. 2. The 
temperature range near t c (nr>) is displayed in the insets; note the different scales. 
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FIG. 4. Phonon correlation function [in 10~ 27 cm 2 s] vs. phonon energy in [meV] for different temperatures below T c (nD) at 
fixed angle 6 = 0. As before, t c (nn) = 0.444, and ( = ka /2 3/2 n = 0.02 was used. 
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FIG. 7. Cluster ordering (precursor) temperature T or d and mean-field (local) transition temperature T c (nu) [in K] as function 
of the defect concentration no; the numerical values of Table I were used here. 
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FIG. 8. Order parameter vs. reduced temperature for systems with different types of defects; the numerical values of Table 
I and no = 10~ 3 were used here. 
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